Transcriptomic landscape of human induced pluripotent stem cell-derived osteogenic differentiation identifies a regulatory role of KLF16.

Osteogenic differentiation is essential for bone development and metabolism, but the underlying gene regulatory networks have not been well investigated. We differentiated mesenchymal stem cells, derived from 20 human induced pluripotent stem cell lines, into preosteoblasts and osteoblasts, and performed systematic RNA-seq analyses of 60 samples for differential gene expression. We noted a highly significant correlation in expression patterns and genomic proximity among transcription factor (TF) and long noncoding RNA (lncRNA) genes. We identified TF-TF regulatory networks, regulatory roles of lncRNAs on their neighboring coding genes for TFs and splicing factors, and differential splicing of TF, lncRNA, and splicing factor genes. TF-TF regulatory and gene co-expression network analyses suggested an inhibitory role of TF KLF16 in osteogenic differentiation. We demonstrate that in vitro overexpression of human KLF16 inhibits osteogenic differentiation and mineralization, and in vivo Klf16+/- mice exhibit increased bone mineral density, trabecular number, and cortical bone area. Thus, our model system highlights the regulatory complexity of osteogenic differentiation and identifies novel osteogenic genes.

The Addis Ababa Lions: Whole-Genome Sequencing of a Rare and Precious Population.

Lions are widely known as charismatic predators that once roamed across the globe, but their populations have been greatly affected by environmental factors and human activities over the last 150 yr. Of particular interest is the Addis Ababa lion population, which has been maintained in captivity at around 20 individuals for over 75 yr, while many wild African lion populations have become extinct. In order to understand the molecular features of this unique population, we conducted a whole-genome sequencing study on 15 Addis Ababa lions and detected 4.5 million distinct genomic variants compared with the reference African lion genome. Using functional annotation, we identified several genes with mutations that potentially impact various traits such as mane color, body size, reproduction, gastrointestinal functions, cardiovascular processes, and sensory perception. These findings offer valuable insights into the genetics of this threatened lion population.

From beer to breadboards: yeast as a force for biological innovation.

The history of yeast Saccharomyces cerevisiae, aka brewer's or baker's yeast, is intertwined with our own. Initially domesticated 8,000 years ago to provide sustenance to our ancestors, for the past 150 years, yeast has served as a model research subject and a platform for technology. In this review, we highlight many ways in which yeast has served to catalyze the fields of functional genomics, genome editing, gene-environment interaction investigation, proteomics, and bioinformatics-emphasizing how yeast has served as a catalyst for innovation. Several possible futures for this model organism in synthetic biology, drug personalization, and multi-omics research are also presented.

Generation of a high yield vaccine backbone for influenza B virus in embryonated chicken eggs.

Influenza B virus (IBV) strains are one of the components of seasonal influenza vaccines in both trivalent and quadrivalent formulations. The vast majority of these vaccines are produced in embryonated chickens' eggs. While optimized backbones for vaccine production in eggs exist and are in use for influenza A viruses, no such backbones exist for IBVs, resulting in unpredictable production yields. To generate an optimal vaccine seed virus backbone, we have compiled a panel of 71 IBV strains from 1940 to present day, representing the known temporal and genetic variability of IBV circulating in humans. This panel contains strains from the B/Victoria/2/87-like lineage, B/Yamagata/16/88-like lineage and the ancestral lineage that preceded their split to provide a diverse set that would help to identify a suitable backbone which can be used in combination with hemagglutinin (HA) and neuraminidase (NA) glycoproteins from any IBV strain to be incorporated into the seasonal vaccine. We have characterized and ranked the growth profiles of the 71 IBV strains and the best performing strains were used for co-infection of eggs, followed by serial passaging to select for high-growth reassortant viruses. After serial passaging, we selected 10 clonal isolates based on their growth profiles assessed by hemagglutination and plaque-forming units. We then generated reverse genetics systems for the three clones that performed best in growth curves. The selected backbones were then used to generate different reassortant viruses with HA/NA combinations from high and low titer yielding wild type IBV. When the growth profiles of the recombinant reassortant viruses were tested, the low titer yielding HA/NA viruses with the selected backbones yielded higher titers similar to those from high titer yielding HA/NA combinations. The use of these IBV backbones with improved replication in eggs might increase yields for the influenza B virus components of seasonal influenza virus vaccines.

The cellular response to drug perturbation is limited: comparison of large-scale chemogenomic fitness signatures.

BACKGROUND: Chemogenomic profiling is a powerful approach for understanding the genome-wide cellular response to small molecules. First developed in Saccharomyces cerevisiae, chemogenomic screens provide direct, unbiased identification of drug target candidates as well as genes required for drug resistance. While many laboratories have performed chemogenomic fitness assays, few have been assessed for reproducibility and accuracy. Here we analyze the two largest independent yeast chemogenomic datasets comprising over 35 million gene-drug interactions and more than 6000 unique chemogenomic profiles; the first from our own academic laboratory (HIPLAB) and the second from the Novartis Institute of Biomedical Research (NIBR). RESULTS: Despite substantial differences in experimental and analytical pipelines, the combined datasets revealed robust chemogenomic response signatures, characterized by gene signatures, enrichment for biological processes and mechanisms of drug action. We previously reported that the cellular response to small molecules is limited and can be described by a network of 45 chemogenomic signatures. In the present study, we show that the majority of these signatures (66%) are also found in the companion dataset, providing further support for their biological relevance as conserved systems-level, small molecule response systems. CONCLUSIONS: Our results demonstrate the robustness of chemogenomic fitness profiling in yeast, while offering guidelines for performing other high-dimensional comparisons including parallel CRISPR screens in mammalian cells.

Human Anti-neuraminidase Antibodies Reduce Airborne Transmission of Clinical Influenza Virus Isolates in the Guinea Pig Model.

The public health burden caused by influenza virus infections is not adequately addressed with existing vaccines and antivirals. Identifying approaches that interfere with human-to-human transmission of influenza viruses remains a pressing need. The importance of neuraminidase (NA) activity for the replication and spread of influenza viruses led us to investigate whether broadly reactive human anti-NA monoclonal antibodies (MAbs) could affect airborne transmission of the virus using the guinea pig model. In that model, infection with recent influenza virus clinical isolates resulted in 100% transmission from inoculated donors to recipients in an airborne transmission setting. Anti-NA MAbs were administered either to the inoculated animals on days 1, 2, and 4 after infection or to the naive contacts on days 2 and 4 after donor infection. Administration of NA-1G01, a broadly cross-reactive anti-NA MAb, to either the donor or recipient reduced transmission of the A/New York City/PV02669/2019 (H1N1) and A/New York City/PV01148/2018 (H3N2) viruses. Administration of 1000-3C05, an anti-N1 MAb, to either the donor or recipient reduced transmission of A/New York City/PV02669/2019 (H1N1) virus but did not reduce transmission of A/New York City/PV01148 (H3N2) virus. Conversely, 229-2C06, an anti-N2 MAb, reduced transmission of A/New York City/PV01148 (H3N2) but did not impact transmission of A/New York City/PV02669/2019 (H1N1) virus. Our work demonstrates that anti-NA MAbs could be further developed into prophylactic or therapeutic agents to prevent influenza virus transmission to control viral spread. IMPORTANCE The burden of influenza remains substantial despite unremitting efforts to reduce the magnitude of seasonal influenza epidemics and prepare for pandemics. Although vaccination remains the mainstay of these efforts, current vaccines are designed to stimulate an immune response against the viral hemagglutinin. Interest in the role immunity against neuraminidase plays in influenza virus infection and transmission has recently surged. Human antibodies that bind broadly to neuraminidases of diverse influenza viruses and protect mice against lethal viral challenge have previously been characterized. Here, we show that three such antibodies inhibit the neuraminidase activity of recent isolates and reduce their airborne transmission in a guinea pig model. In addition to contributing to the accumulating support for incorporating neuraminidase as a vaccine antigen, these findings also demonstrate the potential of direct administration of anti-neuraminidase antibodies to individuals infected with influenza virus and to individuals for postexposure prophylaxis to prevent the spread of influenza virus.

Respiratory Tract Explant Infection Dynamics of Influenza A Virus in California Sea Lions, Northern Elephant Seals, and Rhesus Macaques.

To understand susceptibility of wild California sea lions and Northern elephant seals to influenza A virus (IAV), we developed an ex vivo respiratory explant model and used it to compare infection kinetics for multiple IAV subtypes. We first established the approach using explants from colonized rhesus macaques, a model for human IAV. Trachea, bronchi, and lungs from 11 California sea lions, 2 Northern elephant seals, and 10 rhesus macaques were inoculated within 24 h postmortem with 6 strains representing 4 IAV subtypes. Explants from the 3 species showed similar IAV infection kinetics, with peak viral titers 48 to 72 h post-inoculation that increased by 2 to 4 log10 PFU/explant relative to the inoculum. Immunohistochemistry localized IAV infection to apical epithelial cells. These results demonstrate that respiratory tissue explants from wild marine mammals support IAV infection. In the absence of the ability to perform experimental infections of marine mammals, this ex vivo culture of respiratory tissues mirrors the in vivo environment and serves as a tool to study IAV susceptibility, host range, and tissue tropism. IMPORTANCE Although influenza A virus can infect marine mammals, a dearth of marine mammal cell lines and ethical and logistical challenges prohibiting experimental infections of living marine mammals mean that little is known about IAV infection kinetics in these species. We circumvented these limitations by adapting a respiratory tract explant model first to establish the approach with rhesus macaques and then for use with explants from wild marine mammals euthanized for nonrespiratory medical conditions. We observed that multiple strains representing 4 IAV subtypes infected trachea, bronchi, and lungs of macaques and marine mammals with variable peak titers and kinetics. This ex vivo model can define infection dynamics for IAV in marine mammals. Further, use of explants from animals euthanized for other reasons reduces use of animals in research.

Real-Time Investigation of a Large Nosocomial Influenza A Outbreak Informed by Genomic Epidemiology.

BACKGROUND: Nosocomial respiratory virus outbreaks represent serious public health challenges. Rapid and precise identification of cases and tracing of transmission chains is critical to end outbreaks and to inform prevention measures. METHODS: We combined conventional surveillance with influenza A virus (IAV) genome sequencing to identify and contain a large IAV outbreak in a metropolitan healthcare system. A total of 381 individuals, including 91 inpatients and 290 healthcare workers (HCWs), were included in the investigation. RESULTS: During a 12-day period in early 2019, infection preventionists identified 89 HCWs and 18 inpatients as cases of influenza-like illness (ILI), using an amended definition without the requirement for fever. Sequencing of IAV genomes from available nasopharyngeal specimens identified 66 individuals infected with a nearly identical strain of influenza A H1N1pdm09 (43 HCWs, 17 inpatients, and 6 with unspecified affiliation). All HCWs infected with the outbreak strain had received the seasonal influenza virus vaccination. Characterization of 5 representative outbreak viral isolates did not show antigenic drift. In conjunction with IAV genome sequencing, mining of electronic records pinpointed the origin of the outbreak as a single patient and a few interactions in the emergency department that occurred 1 day prior to the index ILI cluster. CONCLUSIONS: We used precision surveillance to delineate a large nosocomial IAV outbreak, mapping the source of the outbreak to a single patient rather than HCWs as initially assumed based on conventional epidemiology. These findings have important ramifications for more-effective prevention strategies to curb nosocomial respiratory virus outbreaks.

Integrated Transcriptome and Network Analysis Reveals Spatiotemporal Dynamics of Calvarial Suturogenesis.

Craniofacial abnormalities often involve sutures, the growth centers of the skull. To characterize the organization and processes governing their development, we profile the murine frontal suture, a model for sutural growth and fusion, at the tissue- and single-cell level on embryonic days (E)16.5 and E18.5. For the wild-type suture, bulk RNA sequencing (RNA-seq) analysis identifies mesenchyme-, osteogenic front-, and stage-enriched genes and biological processes, as well as alternative splicing events modifying the extracellular matrix. Single-cell RNA-seq analysis distinguishes multiple subpopulations, of which five define a mesenchyme-osteoblast differentiation trajectory and show variation along the anteroposterior axis. Similar analyses of in vivo mouse models of impaired frontal suturogenesis in Saethre-Chotzen and Apert syndromes, Twist1+/- and Fgfr2+/S252W, demonstrate distinct transcriptional changes involving angiogenesis and ribogenesis, respectively. Co-expression network analysis reveals gene expression modules from which we validate key driver genes regulating osteoblast differentiation. Our study provides a global approach to gain insights into suturogenesis.

Influenza A variants with reduced susceptibility to baloxavir isolated from Japanese patients are fit and transmit through respiratory droplets.

Here we report the isolation of the influenza A/H1N1 2009 pandemic (A/H1N1pdm) and A/H3N2 viruses carrying an I38T mutation in the polymerase acidic protein-a mutation that confers reduced susceptibility to baloxavir marboxil-from patients before and after treatment with baloxavir marboxil in Japan. These variants showed replicative abilities and pathogenicity that is similar to those of wild-type isolates in hamsters; they also transmitted efficiently between ferrets by respiratory droplets.

Influenza Virus Polymerase Mutation Stabilizes a Foreign Gene Inserted into the Virus Genome by Enhancing the Transcription/Replication Efficiency of the Modified Segment.

We previously attempted to establish a reporter influenza virus by inserting the gene for the Venus fluorescent protein into the NS segment of influenza A/Puerto Rico/8/34 (PR8, H1N1) virus to yield WT-Venus-PR8. Although the inserted Venus gene was deleted during serial passages of WT-Venus-PR8, we discovered that the PB2-E712D mutation stabilizes the Venus gene. Here, we explored the mechanisms by which Venus gene deletion occurs and how the polymerase mutation stabilizes the Venus gene. Deep sequencing analysis revealed that PB2-E712D does not cause an appreciable change in the mutation rate, suggesting that the stability of the Venus gene is not affected by polymerase fidelity. We found by using quantitative real-time PCR that WT-Venus-PR8 induces high-level interferon beta (IFN-ß) expression. The induction of IFN-ß expression seemed to result from the reduced transcription/replication efficiency of the modified NS segment in WT-Venus-PR8. In contrast, the transcription/replication efficiency of the modified NS segment was enhanced by the PB2-E712D mutation. Loss of the Venus gene in WT-Venus-PR8 appeared to be caused by internal deletions in the NS segment. Moreover, to further our understanding of the Venus stabilization mechanisms, we identified additional amino acid mutations in the virus polymerase complex that stabilize the Venus gene. We found that some of these amino acids are located near the template exit or the product exit of the viral polymerase, suggesting that these amino acids contribute to the stability of the Venus gene by affecting the binding affinity between the polymerase complex and the RNA template and product.IMPORTANCE The reverse genetics method of influenza virus generation has enabled us to generate recombinant viruses bearing modified viral proteins. Recombinant influenza viruses expressing foreign genes have become useful tools in basic research, and such viruses can be utilized as efficient virus vectors or multivalent vaccines. However, the insertion of a foreign gene into the influenza virus genome often impairs virus replication, and the inserted genes are unstable. Elucidation of the mechanisms of foreign gene stabilization will help us to establish useful recombinant influenza viruses.

Circulation of low pathogenic avian influenza (LPAI) viruses in wild birds and poultry in the Netherlands, 2006-2016.

In this study, we explore the circulation of low pathogenic avian influenza (LPAI) viruses in wild birds and poultry in the Netherlands. Surveillance data collected between 2006 and 2016 was used to evaluate subtype diversity, spatiotemporal distribution and genetic relationships between wild bird and poultry viruses. We observed close species-dependent associations among hemagglutinin and neuraminidase subtypes. Not all subtypes detected in wild birds were found in poultry, suggesting transmission to poultry is selective and likely depends on viral factors that determine host range restriction. Subtypes commonly detected in poultry were in wild birds most frequently detected in mallards and geese. Different temporal patterns in virus prevalence were observed between wild bird species. Virus detections in domestic ducks coincided with the prevalence peak in wild ducks, whereas virus detections in other poultry types were made throughout the year. Genetic analysis of the surface genes demonstrated that most poultry viruses were related to locally circulating wild bird viruses, but no direct spatiotemporal link was observed. Results indicate prolonged undetected virus circulation and frequent reassortment events with local and newly introduced viruses within the wild bird population. Increased knowledge on LPAI virus circulation can be used to improve surveillance strategies.

A humanized MDCK cell line for the efficient isolation and propagation of human influenza viruses.

Here, we developed hCK, a Madin-Darby canine kidney (MDCK) cell line that expresses high levels of human influenza virus receptors and low levels of avian virus receptors. hCK cells supported human A/H3N2 influenza virus isolation and growth much more effectively than conventional MDCK or human virus receptor-overexpressing (AX4) cells. A/H3N2 viruses propagated in hCK cells also maintained higher genetic stability than those propagated in MDCK and AX4 cells.

Mandibular dysmorphology due to abnormal embryonic osteogenesis in FGFR2-related craniosynostosis mice.

One diagnostic feature of craniosynostosis syndromes is mandibular dysgenesis. Using three mouse models of Apert, Crouzon and Pfeiffer craniosynostosis syndromes, we investigated how embryonic development of the mandible is affected by fibroblast growth factor receptor 2 (Fgfr2) mutations. Quantitative analysis of skeletal form at birth revealed differences in mandibular morphology between mice carrying Fgfr2 mutations and their littermates that do not carry the mutations. Murine embryos with the mutations associated with Apert syndrome in humans (Fgfr2+/S252W and Fgfr2+/P253R ) showed an increase in the size of the osteogenic anlagen and Meckel's cartilage (MC). Changes in the microarchitecture and mineralization of the developing mandible were visualized using histological staining. The mechanism for mandibular dysgenesis in the Apert Fgfr2+/S252W mouse resulting in the most severe phenotypic effects was further analyzed in detail and found to occur to a lesser degree in the other craniosynostosis mouse models. Laser capture microdissection and RNA-seq analysis revealed transcriptomic changes in mandibular bone at embryonic day 16.5 (E16.5), highlighting increased expression of genes related to osteoclast differentiation and dysregulated genes active in bone mineralization. Increased osteoclastic activity was corroborated by TRAP assay and in situ hybridization of Csf1r and Itgb3 Upregulated expression of Enpp1 and Ank was validated in the mandible of Fgfr2+/S252W embryos, and found to result in elevated inorganic pyrophosphate concentration. Increased proliferation of osteoblasts in the mandible and chondrocytes forming MC was identified in Fgfr2+/S252W embryos at E12.5. These findings provide evidence that FGFR2 gain-of-function mutations differentially affect cartilage formation and intramembranous ossification of dermal bone, contributing to mandibular dysmorphogenesis in craniosynostosis syndromes.This article has an associated First Person interview with the joint first authors of the paper.

Laser Capture Microdissection of Mouse Embryonic Cartilage and Bone for Gene Expression Analysis.

Laser capture microdissection (LCM) is a powerful tool to isolate specific cell types or regions of interest from heterogeneous tissues. The cellular and molecular complexity of skeletal elements increases with development. Tissue heterogeneity, such as at the interface of cartilaginous and osseous elements with each other or with surrounding tissues, is one obstacle to the study of developing cartilage and bone. Our protocol provides a rapid method of tissue processing and isolation of cartilage and bone that yields high quality RNA for gene expression analysis. Fresh frozen tissues of mouse embryos are sectioned and brief cresyl violet staining is used to visualize cartilage and bone with colors distinct from surrounding tissues. Slides are then rapidly dehydrated, and cartilage and bone are isolated subsequently by LCM. The minimization of exposure to aqueous solutions during this process maintains RNA integrity. Mouse Meckel's cartilage and mandibular bone at E16.5 were successfully collected and gene expression analysis showed differential expression of marker genes for osteoblasts, osteocytes, osteoclasts, and chondrocytes. High quality RNA was also isolated from a range of tissues and embryonic ages. This protocol details sample preparation for LCM including cryoembedding, sectioning, staining and dehydrating fresh frozen tissues, and precise isolation of cartilage and bone by LCM resulting in high quality RNA for transcriptomic analysis.

Multiple Incursions and Recurrent Epidemic Fade-Out of H3N2 Canine Influenza A Virus in the United States.

Avian-origin H3N2 canine influenza virus (CIV) transferred to dogs in Asia around 2005, becoming enzootic throughout China and South Korea before reaching the United States in early 2015. To understand the posttransfer evolution and epidemiology of this virus, particularly the cause of recent and ongoing increases in incidence in the United States, we performed an integrated analysis of whole-genome sequence data from 64 newly sequenced viruses and comprehensive surveillance data. This revealed that the circulation of H3N2 CIV within the United States is typified by recurrent epidemic burst-fade-out dynamics driven by multiple introductions of virus from Asia. Although all major viral lineages displayed similar rates of genomic sequence evolution, H3N2 CIV consistently exhibited proportionally more nonsynonymous substitutions per site than those in avian reservoir viruses, which is indicative of a large-scale change in selection pressures. Despite these genotypic differences, we found no evidence of adaptive evolution or increased viral transmission, with epidemiological models indicating a basic reproductive number, R0, of between 1 and 1.5 across nearly all U.S. outbreaks, consistent with maintained but heterogeneous circulation. We propose that CIV's mode of viral circulation may have resulted in evolutionary cul-de-sacs, in which there is little opportunity for the selection of the more transmissible H3N2 CIV phenotypes necessary to enable circulation through a general dog population characterized by widespread contact heterogeneity. CIV must therefore rely on metapopulations of high host density (such as animal shelters and kennels) within the greater dog population and reintroduction from other populations or face complete epidemic extinction.IMPORTANCE The relatively recent appearance of influenza A virus (IAV) epidemics in dogs expands our understanding of IAV host range and ecology, providing useful and relevant models for understanding critical factors involved in viral emergence. Here we integrate viral whole-genome sequence analysis and comprehensive surveillance data to examine the evolution of the emerging avian-origin H3N2 canine influenza virus (CIV), particularly the factors driving ongoing circulation and recent increases in incidence of the virus within the United States. Our results provide a detailed understanding of how H3N2 CIV achieves sustained circulation within the United States despite widespread host contact heterogeneity and recurrent epidemic fade-out. Moreover, our findings suggest that the types and intensities of selection pressures an emerging virus experiences are highly dependent on host population structure and ecology and may inhibit an emerging virus from acquiring sustained epidemic or pandemic circulation.

Avian Influenza Viruses in Wild Birds: Virus Evolution in a Multihost Ecosystem.

Wild ducks and gulls are the major reservoirs for avian influenza A viruses (AIVs). The mechanisms that drive AIV evolution are complex at sites where various duck and gull species from multiple flyways breed, winter, or stage. The Republic of Georgia is located at the intersection of three migratory flyways: the Central Asian flyway, the East Africa/West Asia flyway, and the Black Sea/Mediterranean flyway. For six complete study years (2010 to 2016), we collected AIV samples from various duck and gull species that breed, migrate, and overwinter in Georgia. We found a substantial subtype diversity of viruses that varied in prevalence from year to year. Low-pathogenic AIV (LPAIV) subtypes included H1N1, H2N3, H2N5, H2N7, H3N8, H4N2, H6N2, H7N3, H7N7, H9N1, H9N3, H10N4, H10N7, H11N1, H13N2, H13N6, H13N8, and H16N3, and two highly pathogenic AIVs (HPAIVs) belonging to clade 2.3.4.4, H5N5 and H5N8, were found. Whole-genome phylogenetic trees showed significant host species lineage restriction for nearly all gene segments and significant differences in observed reassortment rates, as defined by quantification of phylogenetic incongruence, and in nucleotide sequence diversity for LPAIVs among different host species. Hemagglutinin clade 2.3.4.4 H5N8 viruses, which circulated in Eurasia during 2014 and 2015, did not reassort, but analysis after their subsequent dissemination during 2016 and 2017 revealed reassortment in all gene segments except NP and NS. Some virus lineages appeared to be unrelated to AIVs in wild bird populations in other regions, with maintenance of local AIVs in Georgia, whereas other lineages showed considerable genetic interrelationships with viruses circulating in other parts of Eurasia and Africa, despite relative undersampling in the area.IMPORTANCE Waterbirds (e.g., gulls and ducks) are natural reservoirs of avian influenza viruses (AIVs) and have been shown to mediate the dispersal of AIVs at intercontinental scales during seasonal migration. The segmented genome of influenza viruses enables viral RNA from different lineages to mix or reassort when two viruses infect the same host. Such reassortant viruses have been identified in most major human influenza pandemics and several poultry outbreaks. Despite their importance, we have only recently begun to understand AIV evolution and reassortment in their natural host reservoirs. This comprehensive study illustrates AIV evolutionary dynamics within a multihost ecosystem at a stopover site where three major migratory flyways intersect. Our analysis of this ecosystem over a 6-year period provides a snapshot of how these viruses are linked to global AIV populations. Understanding the evolution of AIVs in the natural host is imperative to mitigating both the risk of incursion into domestic poultry and the potential risk to mammalian hosts, including humans.